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DOUBLE LOGARITHMIC STABILITY IN THE IDENTIFICATION OF A SCALAR 
POTENTIAL BY A PARTIAL ELLIPTIC DIRICHLET-TO-NEUMANN MAP 

MOURAD CHOULLI, YAVAR KIAN, AND ERIC SOCCORSI 


Abstract. We examine the stability issue in the inverse problem of determining a scalar potential appearing 
in the stationary Schrodinger equation in a bounded domain, from a partial elliptic Dirichlet-to-Neumann 
map. Namely, the Dirichlet data is imposed on the shadowed face of the boundary of the domain and the 
Neumann data is measured on its illuminated face. We establish a log log stability estimate for the Z, 2 -norm 
(resp. the H~ 1 - norm) of fl 1 . for t > 0, and bounded (resp. L 2 ) potentials. 
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1. Introduction 

1.1. Settings and main result. In the present paper Q is a bounded domain of R", n > 3, with C 2 
boundary T. We denote by u(x) the outward unit normal to T computed at x G T. For £ G § ra_1 fixed, we 
introduce the two following subsets of T 

(1.1) r ± (0 = {xGF; ±£ •!/(*)> 0}, 

and denote by F (resp. G ) an open neighborhood of r+(£) (resp. T_(£)) in T. In what follows F + (£) 
(resp. r_(£)) will sometimes be referred to as the ^-shadowed (resp., ^-illuminated) face of T. Next, given 
q G L°°(fl), real-valued, we consider the unbounded self-adioint operator A„ in L 2 (fl), acting on his domain 

D(A q ) = H'(n)nH 2 (n), as 

Aq = —A + q. 

We assume throughout the entire text that 0 is in the resolvent set of A q (i.e. 0 is not in the spectrum of 
A q ) and put 

Q = {q £ L°°(Cl ; R); 0 is not an eigenvalue of A q }. 
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We establish in Section 2 for any q £ Q and g £ H 1 ^ 2 (r) that the boundary value problem (abbreviated 
to BVP in the sequel) 


q / ( -A + <l) u = 0 in ^ 

' " ' \ u = g on r, 

admits a unique transposition solution u £ 77a (U) = {w £ L 2 (U); Aw £ L 2 (fi)} and that the so-called 
Dirichlet-to-Neumann (DN in short) map 


(1.3) A d v u 

is a bounded operator from H~ l ^ 2 {T) into iJ _3 / 2 (r). For qj £ Q, j = 1,2, we denote by Uj the solution to 
m where q-j is substituted for q. Since u = u± — u 2 satisfies 


f (-A + q ± )u = (q 2 - qi)u 2 in Cl, 
| «=0 on r, 


and (q 2 — qi)u 2 £ L 2 (£l), it holds true that u £ D(A qi ). Therefore d u u £ i? 1 / 2 (r) and 
(1.4) k quq2 = A qi - A q2 £ ^(H-^ 2 (T), /f 1/2 (F)), 


hence the operator 

(1-5) A, 1i93 : g £ n S\F) A quq2 (g) ]G , 

is bounded from iJ _1 / 2 (T) fl £'(F), endowed with the norm of Fl~ l / 2 (T), into H 1 ^ 2 (G). We denote by 
||A 91 , 92 || the norm of A gii?2 in ^(i? _1 / 2 (r) PI £'{F), iJ 1 / 2 (G)). 

In the present paper we examine the stability issue in the inverse problem of determining the potential 
q G Q appearing in CuD from the knowledge of A qo ^ q , where qo is an priori known suitable potential of Q. 

Upon denoting by B\ the unit ball of any Banach space X, we may now state the main result of this 
article as follows. 


Theorem 1.1. For any 6 > 0 and t > 0 we may find two constants c > 0 and c > 0, both of them depending 
only on S and t, such that we have 

(!-6) hi-<c(^\\A quq2 \\+ Inc ln||Ag li92 || 

for any q ± ,q 2 £ Qfl AB it » (n) satisfying (q 2 - 9ikn € and 

(!-7) ||«i < cf||A gi , 92 || + Inc ln||A 9lig2 | 

for any q u q 2 £ Q n SB L 2 (fJ ). 


Let us now briefly comment on Theorem 11.1 


Remarks 1.1. (a) We suppose in Theorem II.II that qj , j = 1,2, are real-valued but it is not hard to see that 
the statement can be adapted at the expense of greater technical difficulties, to the case of complex-valued 
potentials. Nevertheless, for the sake of clarity, we shall restrict ourselves to real-valued potentials in the 
remaining part of this text. 

(b) For s > n/2 and e £ (0, s — n/2) we recall from the interpolation theorem jLM l Theorem 12.4, page 73] 
that iP/ 2+e (n) = [H s (£l), H~ 1 ( f2)]g with 0 = (s — (n/2 + e))/(s + 1). Therefore we have 

Iki - 92|U~ ( n) < C'(s)lki - 9211^)1191 - 92||ff-i ( o) < C'(s)^ 1 ""||9i - 92 ||h-i ( o). 

for any qi,q 2 £ Q such that q 2 £ qi + <5-Bfp(n)> and some constant G(s) > 0, depending only on s. From 
this and (hd then follows that 


II 9i 


92 ||L< 


= (fi) < c 


II A, 


Inc 


hi || Ag 



( 1 . 8 ) 
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(c) Fix t £ (0, +oo). Then, arguing as in the derivation of (11.81) . we find two positive constants c and c such 
that the estimate 


9i - 92 Ili 2 (f 2 ) < c ( ||A 91 ,, 2 || + Inc In ||A, 


*-<?!,<?2I 


-1 


holds uniformly in qi,q 2 £ Q obeying 52 G q 1 However in the particular case where dCl is C'^ +1 , 

we point out that the above estimate is weaker than m- 


1.2. State of the art and comments. The celebrated inverse problem of determining q from the knowledge 
of A 9 was first proposed (in a slightly different setting) by Calderon in [C]. The uniqueness issue was treated 
by Sylvester and Uhlmann in [SI Jj and a log-type stability estimate was derived by Alessandrini in [A]. As 
shown by Mandache in [M], this log-type estimate is optimal. 

All the above mentioned results were obtained with the full data, i.e. when measurements are performed 
on the whole boundary T. Taking the Neumann data on r_(£), while the Dirichlet data is imposed on the 
whole boundary T, Bukhgeim and Uhlmann proved in m that partial information of A q still determines 
uniquely the potential. Their result was improved by Kenig, Sjostrand and Uhlmann in [KSIJ] by measuring 
the Dirichlet data on the shadowed face of T and the Neumann data on the illuminated one. Moreover a 
reconstruction result was derived by Nachman and Street in iNSj from the same data as in [KSUI . 

Stability estimates with partial data go back to Heck and Wang’s article IHW| . where the L°°(fl)-norm of 
q is log log stably recovered from A q with partial Neumann data. The same type of estimate was derived in 
|CDR2| . Both papers require that the Dirichlet data be known on the whole boundary. This constraint was 
weakened by Caro, Dos Santos Ferreira and Ruiz in jCDRlj . These authors proved log-log stability of q with 
respect to a partial DN map associated with Dirichlet (resp. Neumann) data measured on a neighborhood 
of U£ g jvF_(£) (resp. U£ g jvF + (£)) where N in a subset of §" _1 . Their result, which is similar to (11.61) . is 
established for the L p -norm, p £ (l,+oo), of bounded and kF A,p -potentials q with A £ (0,1/p]. Therefore 
(11.71) is valid for a wider class of allowable potentials than in ; CDR1 . 

The derivation of Theorem 11.11 relies on complex geometrical optics (CGO in short) solutions to (11.21) and 
the Carleman inequality established by Bukhgeim and Uhlmann in [BU] . These are the two main ingredients 
of the analysis carried out in |HW] , but in contrast to |HWI , we use here this Carleman estimate to construct 
CGO solutions vanishing on a definite part of the boundary T. 

Notice that usual stability estimates in the inverse problem of determining a potential from the full DN 
map are of log type, while they are of log log type for partial data. Indeed, it turns out that the low 
frequencies of the Fourier transform of the potential are bounded uniformly in all directions by the DN map, 
but that this is no longer the case with the partial data. This technical issue can be remedied by using the 
analytic properties of the Fourier transform. The additional log in the stability estimate for the potential 
may thus be seen as the price to pay for recovering this analytic function by its values in a subdomain, which 
is an ill-posed problem. 


1.3. Outline. The paper is organized as follows. In Sect. [2] we introduce the transposition solution associ¬ 
ated with the BVP (11.21) and rigorously define the various DN maps required by the analysis of the inverse 
problem. Sect. [3] is devoted to building CGO solutions that vanish on some part of the boundary T. These 
functions are useful for the proof of Theorem ll.il given in Sect. [4] 


2. Transposition solutions 


In this section we define the transposition solution to the BVP associated with suitable 


data (f,g), 


( 2 . 1 ) 


f (~A + q)u = f in Cl, 
\u = g on T, 


which play a pivotal role in the analysis of the inverse problem carried out in this paper, 
we start by recalling two useful results for the Hilbert space Ha(CI) = {u £ L 2 {Cl)\ A u £ 
with its natural norm 


IMI,HA(n) — (ll M lli 2 (n) + I|Aw|||2(q) j 


1/2 


To this purpose 
L 2 (Cl)} endowed 
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The first result is the following trace theorem , borrowed from |BU1 Lemma 1.1]. 

Lemma 2.1. For j = 0,1, the trace map 

tjU = u £ f^(fi), 

extends to a continuous operator, still denoted by tj, from into ff _ l _1 / 2 (T). Namely, there exists 

Cj > 0, such that the estimate 

llA?' U llff-3-l/ 2 (r) < CjIMIffAlfi)! 

holds for every u £ H^(LT). 

Let us denote by (•, -)j+ 1 / 2 , j = 0,1, the duality pairing between iP +1 / 2 (T) and T), where the 

second argument is conjugated. Then we have the following Generalized Green formula, which can be found 
in [BU| Corollary 1.2], 

Lemma 2.2. Let q be in Then, for any u £ H&(Q) and v £ H 2 (H), we have 

{A-q)uvdx= / u(A-q)vdx+(t 1 u,t 0 v) 3 / 2 -(t 0 u,t 1 v) 1 / 2 . 

Jn Jn 

Let q £ Q. By the usual IL 2 -regularity property for elliptic BVPs (see e.g. 1LM1 Theorem 5.4, page 165]), 
we know that for each / £ L 2 (Ll) and g £ H 3 / 2 (T), there exists a unique solution S q (f,g) £ H 2 (Ll) to (12.11) . 
Moreover the linear operator S q is bounded from L 2 (Ll) x U 3 / 2 (T) into H 2 (LT), i.e. there exists a constant 
C > 0 such that we have 

( 2 - 2 ) ll<S< ? (/,g)||_f/ 2 ( n ) < C (||/||z, 2 (n) + llffll//3/2( r) ) . 

For further reference we put S q p(f) = S q (f, 0) andiS g ,i(g) = <S g (0, g), so we have S q (f, g) = S q fi{f) +5 g ,i(<?) 
from the linearity of S q . 

Next, applying Lemma l2~2l we get for all (/, g) £ L 2 (fl) x H 3 / 2 (T) and v £TL = Hq(£1) D H 2 (Ll) that 


(2.3) 


- / gd v vda(x) + / fvdx = / u(-A + q)vdx, u = S q (f,g). 

Jr Jn Jq 

In view of the left hand side of (|2.3|) we introduce the following continuous anti-linear form on U 

(2.4) £(v) = — j gd v vda{x) + f fvdx , v £ Tl. 

Jr Jn 


In light of (12.21) . the operator L = loS q $ is bounded in L 2 (Vl). Further, with reference to (12.41) . we generalize 
the definition of the anti-linear form t to (f,g) £ TL* x H~ 1 ' 2 { r), upon setting 

£(v) = ~{g,d v v) 1/2 + (f,v), v £ TL, 

where (•, •) denotes the duality pairing between TL and TL* , conjugate linear in its second argument. For any 
h £ L 2 (fl), L{h) = £ (S q ,o(ti)) satisfies 

(2.5) \L(h)\ < ||fl l ||if-i/2( r )||9 v 5q, ) o(/i)||ij-i/2(r) + \\f\\u* ||«5 9 ,o(/i)llw < C (ll.9'lliY- 1 / 2 (r) + ll/llw) ll^llL 2 (np 

according to (& Hence L is a continuous anti-linear form on L 2 (Q). By R.iesz representation theorem, 
there is a unique vector S q (f,g) £ L 2 (Q) such that we have 

-(g,d v S qt0 (h)) 1/2 + (f,S qfi (h)) = [ 5‘(/, ff )Mx, h£L 2 (n). 

Jn 

Bearing in mind that A q is boundedly invertible in L 2 (f2) (since 0 is in the resolvent set of A q ) and that 
S q fi = Af 1 , we obtain upon taking h = A q v in above identity, where v is arbitrary in TL, that 


( 2 . 6 ) 

Moreover 

(2.7) 


entails 


- {g,d u v)i /2 + (f,v) = f S t {f,g)(-A + q)vdx, v£ H. 
Jn 

ll‘5q(/>5)IU 2 (0) < C + II/IIh*) ■ 
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For each (f,g) G T~L* x H 1 / 2 (T), S q (f,g) will be referred to as the transposition solution to the BVP m- 
As a matter of fact we deduce from (12.611 that 

(2.8) (—A + q)S q (f,g) = f in the distributional sense in f l. 

Let us establish now that the transposition solution S q (f,g) coincides with the classical 7J 2 (ff)-solution 
S q (f,g) to (12.11) in the particular case where (f,g) G L 2 {Q) x iL 3 / 2 (T). 

Proposition 2.1. For any ( f,g ) G L 2 (Yl) x H 3 / 2 (T) we have <S*(/,sO = S q {f,g). 

Proof. Put u = S q (f,g). We have u G directly from (12.81) hence tou G H~ l / 2 ( r). Further, given 

ip G H 1 / 2 ( r), we may find such that t\v = tp by the usual extension theorem (see e.g. [LM l Theorem 

8.3, p. 39]). Applying Lemma T2. 21 for such a test function v, we find that 


/ fvdx = / (—A + q)uvdx 


u(-A + q)v dx + (t 0 it, ip) i/ 2 - 


JQ JU JU 

From this and (12.61) then follows that f n fvdx = f n fvdx — (g, ip) 1/2 + (tou, t\v)\/ 2 , which entails 

(g-t 0 u,<p) 1/2 = 0 . 

Since the above identity holds for any ip G U 1/,2 (r) we obtain that tqu = g and hence u = g on 
yields the desired result. 


r. 


This 

□ 


For g G H 1 /, 2 (T) we put S q l (g) = S q (0,g). We have S ql (g) G H&(Q) by (12.81) . with 

l|£g,l(3)llff A (fi) < G (1 + ||g||Loo(o)) ||fl , ||jy-i/2( r ), 

from (12.71) . Hence S q l G <^’(7L _ 1 / 2 (r), H/\(Q)) so we get 

A q = t l0 G ^(H- 1 / 2 ^), H~ 3 / 2 (T)), 

with the aid of Lemma 12.11 

Finally we have «?*(/, 5 ) = <S*(/,0) + <S*(0,ff), by linearity of S*. Put S* fi (f) = S*(f, 0) and S* A (g) = 
S q (0,g). Since S q 0 G 38(H*, and S ql G ^(7L _ 1 / 2 (r), Ha(YI)), it follows from Lemma l2Tl that 

A, = ho Sl^iH- 1 / 2 ^), H~ 3 / 2 (T)). 

Further, as S^^g) = S% A {g) +S qfi (f), with / = -qS^^g), we get that 

A q — Ao T Rq, 

where the operatori? g 

R q : 9 -> hS g ,o ( -qSl^g )) G H^fT), 
is bounded from H~ l / 2 (T) into H l / 2 (Y). 

Remark 2.1. Let us denote by E g the (continuous) DN map 

s, : g G H 3/2 {T) ^ ii S qA (g) G ff 1/2 (r). 

Then, since q is real-valued, we get from the Green formula that 

J gY, q (h)da(x) = J T, q {g)hda(x), g,h £ H 3/2 (T), 

which entails that (E*) ff3/2 ^ r) = E g . On the other hand, we deduce from Lemma HO that E* = A q . 
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3. CGO SOLUTIONS VANISHING ON SOME PART OF THE BOUNDARY 


In this section we build CGO solutions to the Laplace equation appearing in (11.21) . that vanish on a 
prescribed part of the boundary T. This is by means of a suitable Carleman estimate borrowed from [BU| . 
The corresponding result is as follows. 

Proposition 3.1. For S > 0 fixed, let q £ 5B L Let C, r/ £ §" _1 satisfy ( ■ rj = 0 and fix e > 0 so small 
that T1 = T1(C) = {x £ T; C • i/(x) < —e} 0. Then there exists tq = tq(S) > 0, such that for any t > tq 
we may find ip £ L 2 (fiV) obeying ||"^IIi 2 (si) < CV~ 1 / 2 for some constant C > 0 depending only on S, LI and e, 
and such that the function u = e T ^ +lri ' , ' x (l + ip) £ ILa(O) is solution to the BVP 

( (—A + q)u = 0 in O, 

| « = 0 on n. 


Proof. The proof is made of three steps. 

Step 1: A Carleman estimate. For notational simplicity we write T-t instead of T-i-(C), which isdefined in 
ED, and recall from urn that we may find two constants To = to(<5) > 0 and C = C(S) > 0 such the 
estimate 


Ct 2 [ e~ 2 TX <\v\ 2 dx + t [ \C, ■ v{x)\e~ 2 TX <\d v v\ 2 da{x) 

Jn J r + 

< [ e- 2 TX <\(A~q)v\ 2 dx + T f |C-p(x)|e" 2ra: ' c |^u| 2 da(x), 

Jn Jt_ 

holds for all t > tq and v £ PL. Since T±(— CJ) = Tip(C), the above inequality may be equivalently rewritten 
as 


(3.1) 


Ct 2 f e 2rx <\v\ 2 d x + t f |C ■ v{x)\e 2TX <\d v v\ 2 da{x) 

Jn J r_ 

< [ e 2TX ' c |(A - q)v | 2 dx + t f |C • v(x)\e 2TX <\d v v\ 2 da(x), v £ H, r > r 0 . 

Jn J r + 


In view of more compact reformulation of ED we introduce for each real number r the two following scalar 
products: 

(u,v) T = f e 2TX '^uvdx in L 2 (CL), 

Jn 


and 

{(p,ip) T ^,± = [ n(x)e 2 TX< (pipdcr{x), where /i(x) = |C • u(x)|, in L 2 (r ± ). 

J r± 


We denote by L 2 (fil) (resp., L 2 ^(T±)) the space L 2 (Vl) (resp., T 2 (T±)) endowed with the norm || • || T (resp. 
II ■ l|r,/i,±) generated by the scalar product (■, -) T (resp., (•, -) Ti/i ,±). With these notations, the estimate (13.11) 
simply reads 


(3.2) CT 2 \\v\\ 2 r +T\\d u v\\ 2 T ^_ <\\(A-q)v\\ 2 T + T\\d„v\\ 2 T ^ + , v £ Pi, r > r 0 . 

Step 2: Building suitable transposition solutions. Let us identify Lq(LI) = L 2 (Li) (resp., Lq x (T- p) = L 2 (T-i-)) 
with its dual space, so the space dual to L 2 (fil) (resp., L 2 ^( T±))can be identified with L 2 _ t (CL) (resp. 
L 2 _ t i(T±)). Next we consider the operator 

(3.3) P:v£U^{{ A - q)v, d v v ]T+ ) £ L 2 (fil) x L 2 t JT+), 


which is injective by ED- Therefore, for each (f,g) £ L 2 _ r (fil) x L 2 _ t i(T _), the following anti-linear form 


( 3 . 4 ) 


e(wi,w 2 ) = (f,v) + (g, d v v)-, (w 1 ,w 2 )=Pv, v£U, 
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where (•,•) (resp., (-,-)±) denotes the usual scalar product in L 2 (n) (resp., L 2 (r_|-)), is well defined on 
Ran(P). Moreover it holds true for every (wi,w 2 ) = Pv, where v is arbitrary in 'H 1 that 

(3.5) \ e ( w 1 , W 2 )\ < ||/||-T||v||T + ||fl||-r,M-i,-ll^llr,A*,- 



1 /9 

Thus, upon equipping the space L 2 {n) x L 2 /i (T + ) with the norm |||(wi,w; 2 )|||t = (ll^ill? + T ll w; 2 ||?, M , + ) , 

we derive from (13.21) and (13.51) that 

\£{w 1 ,w 2 )\ < C (r- 1 ||/||_ T + T _1/2 ||3ll-T,/i-i,-) (||(A - q)v || 2 + T\\d„v\\ 2 T ^ + ) 1/2 
< C (r- 1 ||/||_ T + T --1/2 ||i?||—7-,/i — i,—) Ill-P^lllr, 

for some constant C > 0 depending only on <5. As a consequence we have 

(3.6) \i{w 1 ,w 2 )\ < C (r' 1 ||/||_ r + T~ 1/2 \\g\\_ T ^-i |||(mi, w 2 )|||t, (wi, w 2 ) £ Ran(P). 

Let us identify the dual space of L 2 (n) x L 2 ^(T+) with L 2 _ T {n) x I?_ t i{T+) endowed with the norm 

11|(zui, 1 ^ 2 )|||—- 7 - = (||wi[|?_ r + r _1 ||w 2 ||?_ ri/i _ li+ ) 1 ^ 2 . Thus, with reference to (13.3D - (13.4D . we deduce from (13.61) 
upon applying Hahn-Banach extension theorem that there exists ( 111 , 112 ) £ L 2 _ T [£l) x L^ T ^_ 1 (r + ) obeying 

(3.7) 111(^2)111-. <c(T- 1 ||/||_ T+ T- 1 / 2 || ff ||_ Ti/i - 1 ,_), 
where C > 0 is the same as in the right hand side of EH), such that we have 


(ui,(A - q)v) + {v 2 ,d v v)+ = (f,v) + {. g,d v v )_, v £ U. 

Bearing in mind that / £ H* and gx r_ — 7’2Xr+ £ P _1 ^ 2 (r), where yr ± is the characteristic function of T± 
in T, the above identity reads 

((-/), v) - (pxr_ -v 2 Xr+,d„v) 1/2 = / v ± (A - q)vdx, v £ H. 

J n 

Therefore, by (12.61) . v\ is the transposition solution to the following IB VP: 


( 3 . 8 ) U-A + q)v 1 = -f inn 

( Vi = gx r_ - V2Xr + on T. 

Step 3: End of the proof. Set p = t(( + ip). The last step of the proof involves picking £ ^(W 1 ; [0,1]) 
such that ip(x) = 1 for x £ Ff_ and supp(tp) DT C T(/ 2 , and considering the transposition solution V\ to the 
BVP (13.81) associated with / = qe p ' x and g = —ipe p ' x . Thus, putting if = e~ p ' x v\ 1 we derive from (13.71) that 

MoHMI —t < c(r- 1 |ke^||_ T +T- 1 /2|| (/3 e^||_^- 1> _) 

< Ct - 1 / 2 ^ t 0 ( 5)- 1 / 2 + ^), 


where |T_ | denotes the (n— l)-dimensional Lebesgue-measure of P_, Further, bearing in mind that p- p = 0, 
it easy to check that u = e p ' x (l + if) £ H^(n) satisfies (—A + q)u = 0 in n. Finally, since v\ = g on T_ by 
(13.81) . we obtain that if = —g) = —1 and consequently that u = 0 on TV This proves the desired result. □ 


4. Proof of Theorem o 

This section contains the proof of Theorem 11.11 which consists of a succession of five lemmas. 
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4.1. A suitable set of Fourier variables. For e > 0 we choose r) = rj(e) £ ( 0 , 7 t/ 2 ) so small that 

(4.1) ((1 — sin #) 2 + 4cos 2 0) 1 ^ < e, 9 £ ( 7 r /2 — 77 , 7 r /2 + 77 ), 

and we define B e as the set of vectors /3 = (/3i ,..., /?„) £ R n with the following spherical coordinates 


(4.2) 


where s £ (0,1), 9\ £ (ir/2 — 77 , 7 r /2 + 77 ), 62 ( 0 , tt/3), G ( 0 , 7 t) for j = 3,..., n — 2, and <f> £ (0, 2-k). Notice 
that B t has positive Lebesgue measure in R” given by: 

pi />7t/2+?7 pir / 3 pir p'K p2-K / n—2 \ 

(4.3) \B e \= / / /•■•// s n ~ l TT sin fc 0 n _!_ fc dad 0 1 ...dfl„_ 2 d^>O. 

J 0 Jn/2—r] JO J 0 ./0 J 0 \/c=l / 

For further reference, we now establish the following result. 

Lemma 4.1. Let T be any orthogonal transformation in R" that maps £ onto e\ = (1, 0,..., 0). Then for 
all e > 0 and k £ T*B e , there exists £ £ S ra_1 satisfying n ■ £ = 0 and |£ — £| < e. 


' Pi 

= s cos 9 \, 


Pi 

= s cos 0j (nU sin 9i^j , j = 2 ,. 

., n — 2 

Pn-1 

= s sin <j) (n^Ii sin 9 k ^j , 


P n 

< 

= s cos £1 (nrii 2 sin 9^ , 



Proof. Let p = (Pi ,..., (5 n ) £ B t be given by (|4.1[) - (|4.2[) . Introduce £ = (£i, £ 2 ,0, 

sin 9\ 


, 0 ), where 


Ci = 


(- 


sin 61 + 


cos 2 0 \ \ 
COS 2 62 J 


1/2 


5 C2 — “ 


COS #1 

cos 6 2 


1 / 2 ' 


(sin 2 0 i + §^) 

Evidently we have £ E and f3 ■ £ = 0. Moreover it follows from m that 

fsln _ 1 12 I cos 2 gi 

ic- ei i 2 = ( * ; c !s^ ^- d 2 + 4 - 2 *<^ 2 - 

1 cos-* 6 2 

Finally, bearing in mind that k = T*(3 and ei = T£, we obtain the desired result upon taking £ = T*£. □ 

4.2. Alessandrini’s identity and consequence. For e > 0 put F e (£) = T \ rf_(£) and G e (£) = T \ 
ri(-£), where ri(±£) is the same as in Proposition 13.11 Since F e (£) = {x £ T; £ ■ v(x) > — e} (resp., 
G e (£) = {a: € P; £ • i/(x) < e}), it holds true that n e >o-F e (£) = { 2 ; £ T; £ • v(x) > 0} = T+(£) (resp., 
n e>0 G e (£) = {x £ P; £ ■ v(x) < 0} = T_(£)). Thus, from the very definitions of F and G, we may choose 
eo = eo(£, F, G) >0 so small that 

(4.4) F 2£ (£) C F and G 2e (£) C G, e G (0,eo]. 

Having said that we turn now to proving the following statement. 

Lemma 4.2. Let T be the same as in Lemma \4-1\ let r £ [to, + 00 ), where r = tq(5) is defined in Proposition 
\3.1[ and let e £ [0,eo). Then, there exists a constant C > 0, depending only on 6 and e, such that the 
estimate 

' <c(e 2dT \\A qi , q2 \\+T- 1 ^ 


fe - qi)e~'*- x dx 


0 ’ 


holds uniformly in k £ rT*B e and r £ (0,2r). Here ||A 9lig2 || is the &(H 1 / 2 (P) n S'(F), H 1 ^ 2 (G))-norm of 
the operator A 9l>g2 defined in (11.31) (11. 5|) . and d = d(Q) = max Ig fj |a;| < 00 . 

Proof. Fix r £ (0,2 t), k £ rT*B t , and let £ be given by Lemma T4.II Pick i £ R n such that t ■ n = £ ■ £ = 0 
(which is possible since n> 3) and \k + £| 2 = \k\ 2 + \t\ 2 = 4r 2 . Set 

Pj = (-iy T (-i K + { : 1)3£ , j = 1,2, 


2 
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>j ) G H<\(Q) be defined in accordance with Proposition 13.11 in such a way that we 

f (—A + qj)itj = 0 in fl, 

l«i=0 onri((2j-3)C), 


< Ct 1/2 . 

Put u = S* l l (toU 2 ) so we have w = u — U 2 = S qit o((qi — 92 )^ 2 ) G 77 and 

(4.7) ii«> = A gi (t 0 u 2 ) - A q2 ( t 0 u 2 ). 

Upon applying the generalized Green formula of Lemma 12.21 with u = u\ and v = w, we derive from m 
that 

(4.8) / (q 2 - qi)uiu 2 dx = (t 0 ui,t 1 w) 1 / 2 . 

J o 

Notice from the second line of (14.51) that the trace toU 2 (resp., ioiti) is supported in F e (£) = r\rf_(£) (resp., 
G e (£) = r \ ri(-O), where n(±£) is defined in Proposition 13.fI Otherwise stated we have F e (f) = {x € 
T; £ • v(x) > —e} (resp., G e (f) = {x € T; C ' u ( x ) < e}) and hence F e (£) C F 2e (£) (resp., G e (f) C G 2e (£)) 
since |£ — £| < e. From this and (14.41) then follows that supp(f 0 "« 2 ) C F (resp., supp(f 0 wi) C G), which 
together with (14.7D - (14.8I) yields 


and let Uj = e Pj ' x (l + 
have 

(4.5) 

and 

(4.6) 


(4-9) 


(<?2 ^ qi)uiu 2 dx 


< II A qi ,q 2 1 | ||iom ||ff-i/2(r) ||ioW2||_H-i/2(r) • 


Moreover, by (14.5l) - (j4.6|) and the very definition of Uj. j = 1,2, we get that 

(4.10) ||io%||ij-i/ 2 (r) < Cj (\\uj\\ L 2 (n) + \\Auj\\ L 2 (n)) < Cj (\\uj\\ L 2 (fy + \\qjUj\\ L 2 (n)) 

< Ce dT (l + t -1 / 2 ). 


Upon possibly substituting max(l, to(<5)) for to(<5) (which does obviously not restrict the generality of the 
above reasoning) we deduce from (I4.9D - (14.101) that 


(4.11) 



qi)uiu 2 dx 


<Ge 2dT ||A 91 , 92 ||. 


Now the desired result follows readily from (14.61) and (14.111) upon taking into account that uiu 2 = e m ' x (l + 

+ □ 

4.3. Bounding the Fourier coefficients. We introduce the function q : K. ra — > R by setting q = (q 2 —qi)xn, 
where denotes the characteristic function of Q in K™. We aim to upper bound the Fourier transform q 
of q , on the unit ball B of K", by means of the following direct generalization of [AEWZl Theorem 4] for 
complex-valued real-analytic functions. 


Theorem 4.1. Assume that the function F : 2B ra — > C is real-analytic and satisfies the condition 

\d a F( K )\ < K g 2 b, aer, 

for some (K , p ) G x (0,1]. Then for any measurable set E C B with positive Lebesgue measure, there 
exist two constants M = M(p, |F|) > 0 and 9 = 0(p , |F|) G (0,1) such that we have 

ll*1U-(B) < MK 1 - 6 ^ jf |F(«)|d^ . 


The result is as follows. 
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Lemma 4.3. For all e £ (0, eo] there exist two constants C = C(Cl,S,e) > 0 and 6 = 8(Q,e) £ (0,1), such 
that we have 

|g(«OI < C'e (1_9)r (e dT \\A qi ,q 2 1 | +r -1/2 ) , k £ rB, re (0,2 r), r e [r 0 ,+cx)). 

Proof. Let r e [to, +oo) and r e (0, 2r) be fixed. By Lemma l4~2l we have 
(4.12) |g(r«)| < C (e 2 dT \\A quq2 1| + r" 1 / 2 ) , « G T*B e . 

In view of (14.31) we apply Theorem 14.11 with E = T*B t and F(k) = q(rn) for n e 2B”. Indeed, in this 
particular case it holds true for any a = (oq,..., a n ) e N n and k £ 2B" that 


d a F(n) = / q(x)(—ir) 




n 

\k= 1 


r dxi... dx n , 


y\ a I 


whence 

(4.13) \d a F( K )\ < \\q\\mn)r'- 'w • ^ ^||a||Li ( n )J - Ul ) 1{a , , 

where we recall that d = d(Q) = max ie fj |x|. Thus, with reference to (I4.12II - (I4.13I) . we obtain that 


W d W < Hall r.irni —rr 


< (iklh 


r\ M ! 
> (rl~l\ 


(4.14) \q{r K )\<Ce^-^ r (e^WAq^qJ + T-^Y, k £ B, 

which immediately yields the result. 


□ 


4.4. Stability inequalities. We turn now to proving the stability inequalities (ll.6l) - (ll.7l) . We start with 
(11.61) . For t > 0 fixed, we assume that q £ F[ t {R n ) and put M = ||a||_f/t(R"). By Parseval inequality, it holds 


true for 

every r > 0 that 



(4.15) 

M L2(n) = 1 llz, 3 (M") 

= 

[ \q(K)\ 2 dk+ f \q{n)\ 2 dk 

J | ac| <r* J | ac| >r* 



< 

f \q(n)\ 2 dk+ f |(1 + H 2 )*|a(«0| 2 dre 

J |«|<r ^ J | ac | >r 



< 

r M 2 

/ m\ 2 dk+ w , 

since J R . 

n l(i + M 2 )*la(«OI 2 dK = ||a||^ t(n) . 

From this and Lemma Pi. 31 then follows that 

(4.16) 

Iklli^n) < 1 

[e dT 

~ \ 2B M 2 

II A 9 i ,92 II + r~ 1/2 J + -%■, r £ (0, 2t), t G [t 0 , + 00 ) 

Let us suppose that 




(4-17) l|A gi ,, 2 ||e(0, 7 o), 

where 70 = v(t 0 ) and v(t) = r _ 1 / 2 e _dT for r £ (0,+oo). Notice that 70 £ (0,1) since r 0 £ [l,+oo). 
Moreover, v being a strictly decreasing function on [ to , +00), there exists a unique r* £ ( to , +00) satisfying 
v >(t*) = ||A gii52 ||. By elementary computation we find that t* = 


T* < 

(2 


(2 | In || Aq 


— In t, 


■) /(2d), 


which entails 


In ||A gii 92 1| , since t* £ ( 1 , + 00 ) and d £ [1, + 00 ). As a consequence the real number t* is greater than 


In || A giig2 1 
(4.18) 


— In 


In ||A gij( j 2 1 


^ /( 2 d), so we get that 

In || A 91j92 I 


T* > 


2 d 



















PARTIAL DN MAP 


11 


upon recalling that In a; < x for all x € (0,+oo). Further, taking r = r* in (14.1611 we obtain for each 
r e (0,2r,) that ||g|| 2 2(n) < C2 2e r n e 2r T~ e + M 2 /r 2t < C'e^ n+2 ^ (2dr*)~ e + M 2 /r 2t , where C' = C2 30 d 0 . 
This and (14.1811 entail 


(4.19) 


1^(0) < C'e^ 


-8 M 2 

ln II A 91,Q2II + —, r e (0,2r*). 


The next step of the derivation involves finding r* € (0,2t*) such that both terms C ,, e ( n + 2 ) r * ln ||A 9li92 j 
and M 2 /r 2t appearing in the right hand side of (14.191) are equal. This can be achieved upon assuming in 

, in such a way that we have 


( M 2 

1 

\ 1 + 2 1-6 

V (n+2)C' 

r ) 

1+2 1-0 

M 2 


IV 


(4.20) [n + 2) ln||A 9l + 

Indeed, we see from (14.181) that the function :r4 (dr) 2t e ( n + 2 ) dr satisfies 


m(2t») > In ||A qii 


q 1,92 I 


0 (n+2)|ln ||A, 


<21 ><22 II | > 


_ 2 t _ 1+2 1 

ln||A 9 i i92 || + (n + 2) ln||A gij92 | 


and hence +2r*) > ( M 2 /C ') ln||A giig2 || according to (14.201) . Thus there exists r* e (0,2r*) such that we 

„ 8 

have id{r*) = (M 2 /C) ln ||A gii92 || . This entails that (2 1 + n + 2)dr if > 2tln(dr*) + (n + 2 )dr* is greater 
that In f(M 2 /C') ln ||A gii92 || J, which combined with the estimate ||g||| 2 (Q) S 2 M 2 /r 2t 1 yields 


(4.21) 


IMIl 2 ( 0 ) < 2 M 


2t + n + 2 

6 


ln 


M 2 \ 

J 


1/8 


In ||A gii q 


Summing up, we obtain (11.61) with c = 2M (( 2 f + n + 2 )/ 6)* and c = (M 2 /C ") 1 / 0 provided ||A 9lj92 || is smaller 
that 7 = minj = o,i 7 j. 

On the other hand, in the particular case where ||A 91j92 || > 7 , we have ||g|| i 2 (f 2 ) < (M/ 7 )||A giig 2 1|. Putting 
this together with (14.211) we end up getting (11.61) upon possibly enlarging c. 

Finally, we obtain (HI) by arguing as in the derivation of (11.61) upon preliminarily substituting the estimate 


(O) = / 1(1 + M 2 ) -1 |g(«0| 2 dK+ f |(l + |+)- 1 |g(«)| 2 d K 

J |«|<r J | ac| >r* 

/ |g(«:)| 2 d k+\[ |g(K)| 2 d K 

J I/d <r J |/d>r 


< 


< 


(K)| 2 dA • 


M 2 


J\k\< r ^ 

for (14.151) . This completes the proof of Theorem 14.II 


5. Application to conductivity problem 


In this section we examine the stability issue in the inverse problem of determining the conductivity 
coefficient er appearing in the system 


(5.1) 


J —div(crVit) = / in fl, 
\ u = g on T, 


from the partial DN map. The strategy is to link this inverse problem to the one studied in the first four 
sections of this paper and then apply Theorem 1 1.1 1 in order to derive a suitable stability estimate for a. 
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Assume that a £ W+AC ={c£ W Al, 00 (ri; R); c(x) > c 0 for some c 0 > 0}. Then for any (/, 3 ) £ L 2 ( fi)x 
if 3 / 2 (r), we know from the standard elliptic theory that EH admits a unique solution S a (f,g) £ H 2 (Gl), 
and that the linear operator 

S* : (f,g) £ L 2 (Q) x H 3/2 (T) H-C(/, 3 ) £ tf 2 (C 

is bounded. In the more general case where (/, g) £ H* xh - 1 / 2 (r), we obtain by arguing in the exact same 
way as in Section [2] that there exists a unique u £ L 2 (Q) obeying 

(5.2) - f udiv(crVTJ)dx = (/, v) - (g, ad u v) 1 / 2 , v GW.. 

J n 

Such a function u will be referred to as the transposition solution of EH) and will be denoted by 5* (/, g). 
Let us introduce the Hilbert space Ciiv(irV)O = {u £ L 2 (f2), div(crVu) £ L 2 ( f2)}, endowed with the 

norm ||w||ij div(o . v) (n) = (|MIz, 2 (q) + ||div(<TVu)|| 2 2 ( n )^ . By a slight modification of the proof of Lemma 

0(e.g. EE]), the trace map 

t°u = cr J 9>|r, u £ j = 0, 1, 

is extended to a linear continuous operator, still denoted by t? , from -ffdiv(o-v)(C into U~ J_ 1 / 2 (r). Thus, 
bearing in mind that S t ( 0 , 3 ) £ 7La(C for 3 £ IL - 1 / 2 (r), we see that the DN map 

A CT : 3 £ if 1 / 2 (r) t?5‘(0, 3 ) £ LT 3 / 2 (r), 


is a bounded operator. 

Assume that er £ W+°°(C = fF 2 ’°°(C 0 Taking into account that 

—div ^ctV(ct _ 1 ^ 2 u)^ = <t 1//2 Au + cr _1 / 2 (Acr 1 / 2 )u^ , v GW, 

we get upon substituting (0, a~ 1 ^ 2 v) for (/, u) in (15.21) . that 

/ o- 1 ^ 2 u(—Au + g CT u)dx =-(cr 1 ^ 2 3 ,v £ "H, 

Jn 

where q a = ct~ 1 / 2 Act 1 / 2 . As a consequence we have A 2 <S* (0,3) = A (0, a Y ^ 2 g), and hence 

o- 1/2 5*(0,ct~ 1/2 3) = A (0,3). 

From this and the identity t\{a x l 2 w) = a~ 1 ^ 2 t1w + \a~ x l 2 [d v a)t qw, which is valid for every w £ H 2 ( fl), 
and generalizes to w £ -ffdiv(o-v) O by duality, we get that 

(5.3) A qa = {d„a)I + cr“ 1/2 A 0 -cr _1/2 . 

Now pick oi ,<72 £ W^’°°(fl) such that o\ = 02 on T and d v cr\ = S„C2 on F. Thus, putting qj = q a . for 
j = 1,2, we deduce from (15.31) that 

(5.4) (A 91 - A 92 ) ( 3 ) = a- 1/2 (A ffl - A CT2 )(ar 1/2 3), 3 £ tf- 1/2 (C n AC- 

Let us next introduce 


A CT1 , CT2 : 3 £ tf-VA) n AC ^ (A ffl - A CT 2 )( 3 )| G £ C/ 2 (G). 


We notice from (15.41) that 


Cl ,q 2 3 _ ^ A(T 1 ,(T2( (7 l ^ 3)1 

for every 3 £ 7L _ 1 / 2 (r) D AC; provided oi = 02 on T and = 
find a constant C > 0, such that we have 


d v d 2 on F D G. In this case, we may 


(5-5) ||A 91i92 || < C'||A CTliCT2 ||, 

where || • |j still denotes the norm of 38 {H~ 1 / 2 {T)r\£'(F ), U 1 / 2 (G)). Here we used the fact that the multiplier 
by cr^ 1 / 2 is an isomorphism of iL ± 1 / 2 (T). 
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Finally, taking into account that <j> = ' 2 — of 2 is solution to the system 

f (-A + = (T 1 2 /2 {q 2 - qi) in 

\ (j) = 0 on r, 

we get ||0 ||l2(q) < C\\q 2 - gi|| ff -i(n) upon taking / = a 1 2 /2 (q 2 - qi) and g = 0 in J2/7]). Thus, applying 
Theorem II.II and recalling (15.51) . we obtain the: 

Corollary 5.1. Let F and G be the same as in section [3 and let S > 0 and <tq > 0. Then for any 
<jj £ SB W 2 ,oa^, j = 1,2, obeying <jj > a o and the condition 

(5.6) a i = <72 on T and d v <j\ — d v u 2 on F fl G, 

we may find a constant C > 0, independent of ay and a 2 , such that we have: 


|ui - 02 ||L 2 (n) < C ||A CTl , CT2 || + 


InC 


In || A 0 



Remark 5.1. It is not clear how to weaken assumption (15.61) in Corollary 15.11 Indeed, to our knowledge, 
the best available result in the mathematical literature (this is a byproduct of |AGI Theorem 2.2, p. 922 
and Theorem 2.4, p. 923]) on the recovery of the conductivity at the boundary, claims for any non empty 
open subset To of T, and for all cr m , m = 1,2, taken as in Corollary 15.11 and satisfying the condition 
supp (cri — a 2 ) |p C r 0 instead m , that 

I- dia 2 IU-^d < C'HA^jVd+A j = 0 , 1 . 

Here A° 1j[T2 denotes the operator g £ H~ X / 2 {T) n £'{T 0 ) ^ (A CTl - A CT2 )(3)| ro £ H^fTo), ||A° 1(T2 || is the 
norm of A ° 1CT2 in D S'(To), iF 1 / 2 (ro)) and C > 0 is a constant that depends neither on a\ nor 

cr 2 . 
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